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Abstract. Observations and iV-body simulations both support a simple 
relation for the disruption time of a cluster as a function of its mass of the 
form: t&\ s = t± x (M / 10 4 Mq) 1 . The scaling factor t$ seems to depend 
strongly on the environment. Predictions and observations show that 
7 ~ 0.64 ±0.06. Assuming that i^is M 0,64 is caused by evaporation and 
^f-) ■ shocking implies a relation between the radius and the mass of a cluster 

CO \ of the form: r^ oc M 0,07 , which has been observed in a few galaxies. The 

suggested relation for the disruption time implies that the lower mass 
end of the cluster initial mass function will be disrupted faster than the 
& • higher mass end, which is needed to evolve a young power law shaped mass 

function into the log-normal mass function of old (globular) clusters. 
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1. Introduction 

This study combines the results of theoretical and observational studies to the 
life time of stellar clusters. It is of vital importance to understand how disrup- 
tion processes influence a cluster population in order to get better understanding 
of the observed luminosity and mass distribution. A key question in studying 
cluster mass functions is whether the power law mass distribution of young clus- 
ters, as observed in interacting/merging and starburst galaxies, will evolve into 
a log-normal mass distribution as observed for old globular cluster populations. 
If so, a preferential depletion of low mass clusters is necessary. The two most 
important dynamical processes influencing the cluster mass function are 1) evap- 
oration, which is due to internal relaxation of the cluster and 2) tidal shocking. 
The first is a result of two body interactions within the cluster which drive the 
low mass stars through interaction with higher mass stars to the edge of the 
cluster until they reach a velocity greater than the escape velocity and leave the 
cluster. Tidal shocks can be caused by passing giant molecular clouds (GMCs) 
or when the cluster travels through a disk, bulge or spiral arm. The stars in the 
cluster gain kinetic energy which causes the cluster to expand and become less 
gravitationally bound. 
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2. Analytical expressions for the disruption time 

The time scales involved for relaxation and shocks depend on the mass and the 
half-mass radius of the cluster: i r h oc M 1 / 2 rj^ 2 and i s h oc Mr t 3 . Boutloukos 
& Lamers (2003) (BL03) and Zepf et al. (1999) have proposed that if a general 
dependence between the disruption time and the mass exists of the form i^is oc 
M 7 , that there should be a relation between the mass and the radius of the 
cluster: oc M 0,07 in order to reduce both time scales to the simple expression 
tih oc ^sh oc M~ 2 / 3 . This shallow relation between radius and mass has been 
observed for clusters in NGC3256 (Zepf et al. 1999) and several spiral galaxies 
(Larsen, these proceedings) 

3. Empirical determination of disruption time 

BL03 propose a scaling law for the disruption time of a star cluster: 

t dis = Ux(M/W 4 M Q y (1) 

with £4 the typical disruption time of a 1O 4 M cluster and 7 a dimensionless 
index. t± and 7 can be derived from the age and mass distributions of a cluster 
sample. BL03 did this for different galaxies and found that £4 is varying per 
galaxy but that 7 is remarkably constant with 7 ~ 0.62 ± 0.06. Disruption 
according to Eq. 1 is needed to explain the observed mass and age distribution 
of clusters in M51 (see also Gieles et al., these proceedings) 

4. Disruption of clusters with iV-body simulations 

Baumgardt & Makino (2003) have studied disruption of clusters of different mass 
and in different orbits in an external tidal field by means of iV-body simulations. 
A clear power-law relation between the disruption time and the initial mass is 
found. The mean predicted slope, corresponding to the 7 value, is 0.66 ± 0.04 
which agrees very well with the empirical value. The scaling value £4 found 
in these simulations depends strongly on the distance to the galactic center 
and ranges between 2 Gyr and 10 Gyr for clusters in a spherical symmetric 
potential. This is much larger then the values found by BL03 (±1 Gyr for the 
solar neighborhood). This can be explained by the different definitions of the 
empirical disruption time and the calculated life time and by the fact that the 
simulations did not include a disk, spiral arms and GMCs which all decrease the 
life time of the cluster in the disk of our galaxy. 
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